Routes to synchrony between asymmetrically interacting oscillator ensembles 
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We report that asymmetrically interacting ensembles of oscillators (AIEOs) follow novel routes 
to synchrony. These routes seem to be a characteristic feature of coupling asymmetry. We show 
that they are unaffected by white noise except that the entrainment frequencies are shifted. The 
probability of occurrence of the routes is determined by phase asymmetry. The identification of 
these phenomena offers new insight into synchrony between oscillator ensembles and suggest new 
ways in which it may be controlled. 
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Coupling asymmetry between oscillators is widespread 
in real physical and biological systems. Examples include 
cardio-respiratory and cardio-<5(EEG) interactions in- 
teractions among activator-inhibitor systems Q , coupled 
circadian oscillators Q , and the interactions between en- 
sembles of oscillators in neuronal dynamics [3, H| . Syn- 
chrony can be sometimes be desirable as in e.g. lasers and 
Josephson- Junction arrays j]], or temporal coding and 
cognition via brain waves [81; but it can also be danger- 
ous, as in epileptic seizures [{|, Parkinson's tremor [lfjj, 
or pedestrians on the Millennium Bridge Thus the 
control of synchronization Q can often be important. An 
understanding how synchrony arises is of course an essen- 
tial prerequisite for the development of control schemes. 

In this Letter, we report a phenomenon that occurs 
in asymmetrically interacting ensembles of oscillators 
(AIEOs), where it provides novel routes to synchrony. 
These routes are the characteristic of AIEOs and they 
cannot occur in systems where the interactions are sym- 
metrical. Our results yield new insights into how syn- 
chrony arises and offer possible ways of controlling syn- 
chrony between real AIEOs. 

The phase dynamical equations of a system of two 
AIEOs can be written as \\% 
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cies u>l ' are assumed to be Lorentzianly distributed 
as g(cj( 1,2 ') = -(7 2 + {to — a^ 1 ' 2 )) 2 ) -1 with central fre- 
quencies u) 1,2 , and 7 is the half-width at half-maximum. 
With this characterization, we show that an increase of 
the coupling strength between two ensembles that are 
synchronized separately does not immediately result in 
their mutual phase-locking. Rather, phase-locking oc- 
curs through either one of two different routes: in Route- 

1 the oscillators in the two ensembles combine and form 
clusters; in Route-II one of the ensembles desynchronizes 
while the other remains synchronized. Further, there also 
exists the possibility that phase-locking between the en- 
sembles cannot occur at all. 

The model can conveniently be expressed in 
terms of (complex-valued, mean-field) order parameters 
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. Here ^^(t) are the 
average phases of the oscillators in the respective en- 
sembles and r^ l,2 '(t) provide measures of the coherence 
of each oscillator ensemble which varies from to 1. 
When r^ 1 ' 2 ) rj 1 the corresponding ensemble is synchro- 
nized in phase (microscopic synchronization) and when 
S'lp = ip^ — ~ constant the ensembles are mutu- 
ally locked in phase (macroscopic synchronization) . With 
these definitions, Eq. (TT|) becomes 
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The interactions are characterized by coupling param- 
eters A^- 1 ' 2 ^ and B to quantify respectively the inter- 
actions within, and between, the ensembles. The fact 
that AW ^ A^ (c-asymmetry hereafter) implies that 
the oscillators in the ensembles are asymmetrically cou- 

(1 2) 

pled. 9\ ' are the phases of the ith oscillator in each 
ensemble and N^ 1 ' 2 ^ refer to the ensemble sizes; we take 
N (i) _ N (2) _ N ^ with TV — > oo. Phase asymmetry 

(p-asymmetry hereafter) is introduced by phase shifts 



where for simplicity we have considered the particular 
case q^ 1 ' 2 ' 3 ) = a. In the limit N — > oo, a density func- 
tion can be defined as p^ 1 ' 2 ^ {9,t,uj)dtod9 which describes 
the number of oscillators with natural frequencies within 
[cj, lo + cLu] and with phases within [6, 9 + dO] at time 
t. For fixed u the distribution p( 1,2 ) (0,t,u>) obeys the 
evolution equation dpW/dt = ~d(p^9^^)/d9. The 
function p( 1,2 > (9, t, lj) is real and 27r periodic in 9, so it 
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FIG. 1: (B—j) bifurcation diagram for a = 0, Alu — 1. The 
different synchronization regimes are described in the text. 
The boundary between regimes NS and S/D represents 7 C +. 
The *s represent the numerical bifurcation boundaries; (left) 
^(i) = A {2) = 1, (right) A w = 1.2,yl (2) = 1. Insets show the 
frequency distributions (also obtained numerically) for the in- 
dicated regions; their ordinate axes represent oscillator counts 
in thousands. Note that the occurrence of perfect synchro- 
nization with 2000 (left) and 1000 (right) oscillator groups 
will not occur throughout all of each indicated region. 
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where c.c is the complex conjugate of the preceding term 
and r)(6,t,uj) denotes the 2nd and higher harmonics. 
Substituting p( 1,2 ) (9, t, u>) into the evolution equation, we 
get the following linearized equation for p\ as 



Pi 



(1,2) 



■ (1,2) 

nop\ 



where the Fourier components for \l\ > 2 are neglected 
since I = ±1 are the only nontrivial unstable modes and 
po = 1/27T is the trivial solution corresponding to inco- 
herence. (■) represents the average over the frequencies 
w (i,2) weighted by the Lorentzian distribution g(a/ 1,2 )). 
The eigenvalues obtained from the characteristic equa- 
tion of ((3]) are 
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where k = A«+A( 2 ),£ = (\A 2 + B 2 ),A= (A^-A^), 
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p = ^cos(2a) + AAuJsma — Aw 2 , and q — £sin(2a) — 
AAu> cos a. As a signature of synchronization, we take 
the condition Re(X±) > for analytic treatment. For 
the numerical experiment, we set A 1 ' 2 ) > 0.7 for mi- 
croscopic synchronization in the corresponding ensembles 
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FIG. 2: (B-a) bifurcation diagram for (a) A {1) = 1.2, A (2) = 
1, (b) A {1) = 1.8, A^ = 1.4 and Atu = 1,7 = 0.5. Note that 
the synchronization regimes SI and Dl are greatly reduced in 
the presence of p-asymmetry. The *s denote the numerical 
bifurcation boundary between the synchronized and incoher- 
ent states. The discrepancy between numerical and analytical 
boundaries occurs due the p-asymmetry (unlike Fig. [TJ that 
affects region S2 and hence the thresholds for r' 1 ' and r' 2 '. 



and a constant Sip for macroscopic synchronization as the 
conditions. We take TV = 1000 in each ensemble and the 
equations are solved using a fourth-order Runge-Kutta 
routine. The initial phases of the oscillators are assumed 
to be equally distributed in [0, 2w). 

One might intuitively anticipate the possibility of four 
distinct dynamical regimes: no synchronization (NS); 
global synchronization, in which the oscillators of both 
ensembles are entrained to the same frequency (SI); syn- 
chronization within one ensemble but not the other (S2); 
synchronization within both ensembles, separately and 
independently, with two entrainment frequencies (D2). 
In what follows, we show that this is indeed the case 
and, furthermore, that there is a global regime in which 
the two ensembles behave as one, but oscillators within 
each ensemble are entrained at either one of two distinct 
entrainment frequencies (Dl). Regions S2 and Dl cannot 
occur when c-asymmetry is absent 13( (see Fig. [TJ. 

For the case a = 0, when p > 0, in the region SI, the 
incoherent (steady) state becomes unstable via a single 
Hopf bifurcation and the ensembles entrain to a single fre- 
quency f2+. With further decrease of 7 below the 7 C _ line 
in the region Dl, a new entrainment frequency emerges 
through a second Hopf bifurcation. In this region, the os- 
cillators from the two ensembles combine and form two 
clusters (macroscopic clustering) oscillating with two fre- 
quencies, Q± = -Im(A±) = ±{l/2){p 2 + q 2 ]i sin(±()+oj. 
7 C ± are obtained by imposing the condition i?e(A±) = 0. 
Thus in this region the order parameters either fluctuate 
in a quasi-periodic manner or have complicated dynam- 
ics (see Fig. [3]). This is because each ensemble has two 
clusters oscillating with different frequencies (see Figs. 
[3ja) and [4]) . The presence of two entrainment fre- 
quencies can be seen by looking at the frequencies into 
which all the individual oscillators are grouped as shown 
in Fig. [5] The macroscopic clustering that occurs in 
this case is quite different from formation of clusters in a 
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FIG. 4: Histograms showing virtual desynchronization - Left: 
First ensemble, Right: Second ensemble, for the correspond- 
ing parameter values in Fig. [3] (a) . Note that the first fre- 
quency component in both the ensembles have two indistin- 
guishably different subcomponents which may be considered 
as one. This occurs due to the discrepancy between numerics 
and analytics. 



FIG. 3: Numerical plot of coherence parameters (grey), 
r' 2 ' (black) and phase difference Sip as a function of time. 
Here (a) B = 1, a = 0.23 and (b) B = 1, a = OAT. (a) and 
(b) correspond to regions DI and D2 of Fig. [2Jb) for the same 
values of parameters as traveling along the line I — II. Note 
that the order parameters display no synchronization in the 
Dl region. 



single ensemble IToT j - here the oscillators in two dif- 
ferent ensembles combine and form clusters. The occur- 
rence of this phenomenon provides new insight into the 
control of synchrony in realistic situations where there 
is asymmetry, like neural networks where neurons from 
one ensemble (e.g. cortex) tend to synchronize with those 
in the other ensemble (e.g. thalamus) thus giving rise to 
creating desirable or undesirable (e.g. epileptic seizures) 
effects. In the absence of coupling asymmetry, these phe- 
nomena do not exist. Traversing the line I— II of Fig. Q] 
demonstrates Route-I to phase-locking of the ensembles. 
As we increase B, the oscillators in the ensembles pass 
from the dynamical state of microscopic synchronization 
(D2) through macroscopic clustering to macroscopic syn- 
chronization or phase-locking of the ensembles. Note 
that when A^> = or Auj = only one entrainment 
frequency exists below 7 C _ and therefore this sequence 
does not occur (due to the absence of region Dl). While 
traversing the line III-IV of Fig. [1] the ensembles adopt 
Route-II to phase-locking in which it is the macroscopic 
clustering that does not occur. 

When p < 0, corresponding to regions S2 and D2 in 
Fig. [TJ the dynamics is the same as above except that 
the values of critical T and synchronization frequencies 
differ, as can be calculated from Eq. ((4|). In region S2, 
microscopic synchronization can occur in either one of 
the ensembles, depending upon whether A^ or A^ is 
greater; in Fig. [TJ since > A^ 2 \ synchronization 

occurs in the first ensemble with the second ensemble re- 
maining incoherent. Note that, on increasing B while in 
region S2, the condition p < is violated and the en- 
sembles enter into the phase- locked region SI. In region 



D2, the ensembles synchronize separately to two lock- 
ing frequencies (unlike in region Dl where the ensembles 
combine). The corresponding (B — 7) bifurcation dia- 
gram for the case A(i)=AW is plotted in Fig. Q] (left) to 
show the difference between these two cases. The region 
D represents microscopic synchronization which occurs 
through a degenerate Hopf bifurcation (similar to D2) 
and S represents macroscopic synchronization through a 
single Hopf bifurcation (similar to SI). The regions S2 
and Dl cannot arise for this case. Hence the ensembles 
cannot follow either Route-I or Route-II to synchrony but 
can only pass from the dynamical states of microscopic 
synchrony to macroscopic synchrony. 

For the case a ^ 0, regions Dl and SI shrink as 
a increases, whereas S2 expands, as shown in Fig. [2] 
This means that p-asymmetry reduces the probability of 
macroscopic synchronization (reduced SI and Dl regions 
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FIG. 5: Numerical plot of coherence parameters r ( ' (grey), 
r' 2 ' (black) and phase difference 5tp as a function of time, 
(a) B = 0.7, a = 0.2 (b) B = 0.7, a = tt/4 (near a c ) (c) 
B — 0.7, a = 1.2 (a > a c ). (a), (b) and (c) correspond to 
regions D2, near NS/S2 boundary and NS of Fig. [2ja) for the 
same values of parameters. 
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in Fig. [2]) and mostly allows only microscopic synchro- 
nization of one or both of the ensembles. For a given set 
of parameters, there exists a value of a below which the 
condition p > is satisfied and above which p < is satis- 
fied. As a result, when a > ctj the macroscopic synchrony 
breaks and the system enters the microscopically syn- 
chronized state. Thus as one travels from SI (Dl) to S2 
(D2) across <x, the combined synchrony with single (dou- 
ble) frequency breaks between the ensembles and inde- 
pendent synchronization with single (double) frequency 
regime appears. Region S2, unlike in Fig. [TJ embraces 
two states (i) synchronization in ensemble 1 with ensem- 
ble 2 incoherent and (ii) synchronization in ensemble 2 
with ensemble 1 incoherent, but does distinguish between 
them. Further, there is a critical value of a = a c above 
which the collective oscillations disappear and the inco- 
herent state becomes stabilized (see Figs. [2] and [SJ . 

Real physical systems are of course subject to noise 
(random fluctuations, of either internal or external ori- 
gin), so we now consider how the above analysis will be 
modified as a result. Adding i)f ' 2) to the RHS of Eq. Q, 

(12). 

where 77) ' are independent Gaussian white noises with 
fa (1,2) (t)) = Oand (#' 2) (t)*7f 2) '(*)) - 2K^5(t-t')S l3 
and KM are the noise intensities, the eigenvalues of the 
linearized equation then take the form 
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where K = ( g(1) + g<2) ) ; AK = K^' - A y "> p 
£ cos(2a) + A[Aujsma + AA'cos(a)] - Auj 2 + AK 2 and 
q = £ sin(2a) - AAoj cos a - AAK sin a + 2AujAK. For 
simplicity, if we consider = K^ = K, one can then 
replace 7 in Eq. (TJ| by 7 + -ft'. Thus from Eq. flS} it is 
evident that the dynamics and hence the new route to 
synchrony are unaffected by the presence of white noise. 
Consequently, for increasing noise intensity, the incoher- 
ent state becomes unstable for larger values of the criti- 
cal parameters. When K^' ^ K&> the entrainment fre- 
quency depends on AK (shown in Fig. [6|) . On increasing 
AK (K^ > K^), the microscopic synchronization of 
ensemble 1 is destroyed. 

In summary, we have found new routes to synchrony 
between two AIEOs, characterized by asymmetry in the 
coupling strength. We have also established that the ef- 
fect of white noise is simply to alter the critical values of 
the parameters that control bifurcation, and to change 
the entrainment frequencies, but without otherwise af- 
fecting this route to synchrony. Our main focus has been 
on phase oscillator ensembles, but similar phenomena can 
also be identified in ensembles of other kinds of oscilla- 
tors. 
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FIG. 6: The shift in the entrainment frequencies 0± with 
AK. Here A {1) = A {2 \ B = 0.5, a = 0, Q = and Aw = 1. 
N represents oscillator count in thousands. 
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in real physical and biological systems. Examples include 
cardio-respiratory and cardio-<5(EEG) interactions in- 
teractions among activator-inhibitor systems Q , coupled 
circadian oscillators Q , and the interactions between en- 
sembles of oscillators in neuronal dynamics [3, H| . Syn- 
chrony can be sometimes be desirable as in e.g. lasers and 
Josephson- Junction arrays j]], or temporal coding and 
cognition via brain waves [81; but it can also be danger- 
ous, as in epileptic seizures [{|, Parkinson's tremor [lfjj, 
or pedestrians on the Millennium Bridge Thus the 
control of synchronization Q can often be important. An 
understanding how synchrony arises is of course an essen- 
tial prerequisite for the development of control schemes. 

In this Letter, we report a phenomenon that occurs 
in asymmetrically interacting ensembles of oscillators 
(AIEOs), where it provides novel routes to synchrony. 
These routes are the characteristic of AIEOs and they 
cannot occur in systems where the interactions are sym- 
metrical. Our results yield new insights into how syn- 
chrony arises and offer possible ways of controlling syn- 
chrony between real AIEOs. 

The phase dynamical equations of a system of two 
AIEOs can be written as \\% 
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as g(cj( 1,2 ') = -(7 2 + {to — a^ 1 ' 2 )) 2 ) -1 with central fre- 
quencies u) 1,2 , and 7 is the half-width at half-maximum. 
With this characterization, we show that an increase of 
the coupling strength between two ensembles that are 
synchronized separately does not immediately result in 
their mutual phase-locking. Rather, phase-locking oc- 
curs through either one of two different routes: in Route- 

1 the oscillators in the two ensembles combine and form 
clusters; in Route-II one of the ensembles desynchronizes 
while the other remains synchronized. Further, there also 
exists the possibility that phase-locking between the en- 
sembles cannot occur at all. 

The model can conveniently be expressed in 
terms of (complex-valued, mean-field) order parameters 
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of each oscillator ensemble which varies from to 1. 
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S'lp = ip^ — ~ constant the ensembles are mutu- 
ally locked in phase (macroscopic synchronization) . With 
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The interactions are characterized by coupling param- 
eters A^- 1 ' 2 ^ and B to quantify respectively the inter- 
actions within, and between, the ensembles. The fact 
that AW ^ A^ (c-asymmetry hereafter) implies that 
the oscillators in the ensembles are asymmetrically cou- 
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ensemble and N^ 1 ' 2 ^ refer to the ensemble sizes; we take 
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where for simplicity we have considered the particular 
case q^ 1 ' 2 ' 3 ) = a. In the limit N — > oo, a density func- 
tion can be defined as p^ 1 ' 2 ^ {9,t,uj)dtod9 which describes 
the number of oscillators with natural frequencies within 
[cj, lo + cLu] and with phases within [6, 9 + dO] at time 
t. For fixed u the distribution p( 1,2 ) (0,t,u>) obeys the 
evolution equation dpW/dt = ~d(p^9^^)/d9. The 
function p( 1,2 > (9, t, lj) is real and 27r periodic in 9, so it 
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FIG. 1: (B—j) bifurcation diagram for a = 0, Alu — 1. The 
different synchronization regimes are described in the text. 
The boundary between regimes NS and S/D represents 7 C +. 
The *s represent the numerical bifurcation boundaries; (left) 
^(i) = A {2) = 1, (right) A w = 1.2,yl (2) = 1. Insets show the 
frequency distributions (also obtained numerically) for the in- 
dicated regions; their ordinate axes represent oscillator counts 
in thousands. Note that the occurrence of perfect synchro- 
nization with 2000 (left) and 1000 (right) oscillator groups 
will not occur throughout all of each indicated region. 



can be expressed as a Fourier series in 9 
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+ c.c + r)(6,t,cj), 



where c.c is the complex conjugate of the preceding term 
and r)(6,t,uj) denotes the 2nd and higher harmonics. 
Substituting p( 1,2 ) (9, t, u>) into the evolution equation, we 
get the following linearized equation for p\ as 



Pi 



(1,2) 



■ (1,2) 

nop\ 



where the Fourier components for \l\ > 2 are neglected 
since I = ±1 are the only nontrivial unstable modes and 
po = 1/27T is the trivial solution corresponding to inco- 
herence. (■) represents the average over the frequencies 
w (i,2) weighted by the Lorentzian distribution g(a/ 1,2 )). 
The eigenvalues obtained from the characteristic equa- 
tion of ((3]) are 



±±( P 2 + q 2 )ie^ 



qc 



A± = 4 



;(p 2 + q 2 ) 



4 e 2 



qc 



p>0 

- iuj, 
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(4) 



where k = A«+A( 2 ),£ = (\A 2 + B 2 ),A= (A^-A^), 



(1) 



;(2) 



)/2, 



p = ^cos(2a) + AAuJsma — Aw 2 , and q — £sin(2a) — 
AAu> cos a. As a signature of synchronization, we take 
the condition Re(X±) > for analytic treatment. For 
the numerical experiment, we set A 1 ' 2 ) > 0.7 for mi- 
croscopic synchronization in the corresponding ensembles 
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FIG. 2: (B-a) bifurcation diagram for (a) A {1) = 1.2, A (2) = 
1, (b) A {1) = 1.8, A^ = 1.4 and Atu = 1,7 = 0.5. Note that 
the synchronization regimes SI and Dl are greatly reduced in 
the presence of p-asymmetry. The *s denote the numerical 
bifurcation boundary between the synchronized and incoher- 
ent states. The discrepancy between numerical and analytical 
boundaries occurs due the p-asymmetry (unlike Fig. [TJ that 
affects region S2 and hence the thresholds for r' 1 ' and r' 2 '. 



and a constant Sip for macroscopic synchronization as the 
conditions. We take TV = 1000 in each ensemble and the 
equations are solved using a fourth-order Runge-Kutta 
routine. The initial phases of the oscillators are assumed 
to be equally distributed in [0, 2w). 

One might intuitively anticipate the possibility of four 
distinct dynamical regimes: no synchronization (NS); 
global synchronization, in which the oscillators of both 
ensembles are entrained to the same frequency (SI); syn- 
chronization within one ensemble but not the other (S2); 
synchronization within both ensembles, separately and 
independently, with two entrainment frequencies (D2). 
In what follows, we show that this is indeed the case 
and, furthermore, that there is a global regime in which 
the two ensembles behave as one, but oscillators within 
each ensemble are entrained at either one of two distinct 
entrainment frequencies (Dl). Regions S2 and Dl cannot 
occur when c-asymmetry is absent 13( (see Fig. [TJ. 

For the case a = 0, when p > 0, in the region SI, the 
incoherent (steady) state becomes unstable via a single 
Hopf bifurcation and the ensembles entrain to a single fre- 
quency f2+. With further decrease of 7 below the 7 C _ line 
in the region Dl, a new entrainment frequency emerges 
through a second Hopf bifurcation. In this region, the os- 
cillators from the two ensembles combine and form two 
clusters (macroscopic clustering) oscillating with two fre- 
quencies, Q± = -Im(A±) = ±{l/2){p 2 + q 2 ]i sin(±()+oj. 
7 C ± are obtained by imposing the condition i?e(A±) = 0. 
Thus in this region the order parameters either fluctuate 
in a quasi-periodic manner or have complicated dynam- 
ics (see Fig. [3]). This is because each ensemble has two 
clusters oscillating with different frequencies (see Figs. 
[3ja) and [4]) . The presence of two entrainment fre- 
quencies can be seen by looking at the frequencies into 
which all the individual oscillators are grouped as shown 
in Fig. [5] The macroscopic clustering that occurs in 
this case is quite different from formation of clusters in a 
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FIG. 4: Histograms showing virtual desynchronization - Left: 
First ensemble, Right: Second ensemble, for the correspond- 
ing parameter values in Fig. [3] (a) . Note that the first fre- 
quency component in both the ensembles have two indistin- 
guishably different subcomponents which may be considered 
as one. This occurs due to the discrepancy between numerics 
and analytics. 



FIG. 3: Numerical plot of coherence parameters (grey), 
r' 2 ' (black) and phase difference Sip as a function of time. 
Here (a) B = 1, a = 0.23 and (b) B = 1, a = OAT. (a) and 
(b) correspond to regions DI and D2 of Fig. [2Jb) for the same 
values of parameters as traveling along the line I — II. Note 
that the order parameters display no synchronization in the 
Dl region. 



single ensemble IToT j - here the oscillators in two dif- 
ferent ensembles combine and form clusters. The occur- 
rence of this phenomenon provides new insight into the 
control of synchrony in realistic situations where there 
is asymmetry, like neural networks where neurons from 
one ensemble (e.g. cortex) tend to synchronize with those 
in the other ensemble (e.g. thalamus) thus giving rise to 
creating desirable or undesirable (e.g. epileptic seizures) 
effects. In the absence of coupling asymmetry, these phe- 
nomena do not exist. Traversing the line I— II of Fig. Q] 
demonstrates Route-I to phase-locking of the ensembles. 
As we increase B, the oscillators in the ensembles pass 
from the dynamical state of microscopic synchronization 
(D2) through macroscopic clustering to macroscopic syn- 
chronization or phase-locking of the ensembles. Note 
that when A^> = or Auj = only one entrainment 
frequency exists below 7 C _ and therefore this sequence 
does not occur (due to the absence of region Dl). While 
traversing the line III-IV of Fig. [1] the ensembles adopt 
Route-II to phase-locking in which it is the macroscopic 
clustering that does not occur. 

When p < 0, corresponding to regions S2 and D2 in 
Fig. [TJ the dynamics is the same as above except that 
the values of critical T and synchronization frequencies 
differ, as can be calculated from Eq. ((4|). In region S2, 
microscopic synchronization can occur in either one of 
the ensembles, depending upon whether A^ or A^ is 
greater; in Fig. [TJ since > A^ 2 \ synchronization 

occurs in the first ensemble with the second ensemble re- 
maining incoherent. Note that, on increasing B while in 
region S2, the condition p < is violated and the en- 
sembles enter into the phase- locked region SI. In region 



D2, the ensembles synchronize separately to two lock- 
ing frequencies (unlike in region Dl where the ensembles 
combine). The corresponding (B — 7) bifurcation dia- 
gram for the case A(i)=AW is plotted in Fig. Q] (left) to 
show the difference between these two cases. The region 
D represents microscopic synchronization which occurs 
through a degenerate Hopf bifurcation (similar to D2) 
and S represents macroscopic synchronization through a 
single Hopf bifurcation (similar to SI). The regions S2 
and Dl cannot arise for this case. Hence the ensembles 
cannot follow either Route-I or Route-II to synchrony but 
can only pass from the dynamical states of microscopic 
synchrony to macroscopic synchrony. 

For the case a ^ 0, regions Dl and SI shrink as 
a increases, whereas S2 expands, as shown in Fig. [2] 
This means that p-asymmetry reduces the probability of 
macroscopic synchronization (reduced SI and Dl regions 
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FIG. 5: Numerical plot of coherence parameters r ( ' (grey), 
r' 2 ' (black) and phase difference 5tp as a function of time, 
(a) B = 0.7, a = 0.2 (b) B = 0.7, a = tt/4 (near a c ) (c) 
B — 0.7, a = 1.2 (a > a c ). (a), (b) and (c) correspond to 
regions D2, near NS/S2 boundary and NS of Fig. [2ja) for the 
same values of parameters. 
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in Fig. [2]) and mostly allows only microscopic synchro- 
nization of one or both of the ensembles. For a given set 
of parameters, there exists a value of a below which the 
condition p > is satisfied and above which p < is satis- 
fied. As a result, when a > ctj the macroscopic synchrony 
breaks and the system enters the microscopically syn- 
chronized state. Thus as one travels from SI (Dl) to S2 
(D2) across <x, the combined synchrony with single (dou- 
ble) frequency breaks between the ensembles and inde- 
pendent synchronization with single (double) frequency 
regime appears. Region S2, unlike in Fig. [TJ embraces 
two states (i) synchronization in ensemble 1 with ensem- 
ble 2 incoherent and (ii) synchronization in ensemble 2 
with ensemble 1 incoherent, but does distinguish between 
them. Further, there is a critical value of a = a c above 
which the collective oscillations disappear and the inco- 
herent state becomes stabilized (see Figs. [2] and [SJ . 

Real physical systems are of course subject to noise 
(random fluctuations, of either internal or external ori- 
gin), so we now consider how the above analysis will be 
modified as a result. Adding i)f ' 2) to the RHS of Eq. Q, 

(12). 

where 77) ' are independent Gaussian white noises with 
fa (1,2) (t)) = Oand (#' 2) (t)*7f 2) '(*)) - 2K^5(t-t')S l3 
and KM are the noise intensities, the eigenvalues of the 
linearized equation then take the form 



A K=0 



A K=0.5 



-K -l+je ia ±Up 2 + q 2 



p>0 



-K-J+ f e ia ± §(p 2 -I- q 2 )ie 



(5) 



1 2 C — it 
p < 



where K = ( g(1) + g<2) ) ; AK = K^' - A y "> p 
£ cos(2a) + A[Aujsma + AA'cos(a)] - Auj 2 + AK 2 and 
q = £ sin(2a) - AAoj cos a - AAK sin a + 2AujAK. For 
simplicity, if we consider = K^ = K, one can then 
replace 7 in Eq. (TJ| by 7 + -ft'. Thus from Eq. flS} it is 
evident that the dynamics and hence the new route to 
synchrony are unaffected by the presence of white noise. 
Consequently, for increasing noise intensity, the incoher- 
ent state becomes unstable for larger values of the criti- 
cal parameters. When K^' ^ K&> the entrainment fre- 
quency depends on AK (shown in Fig. [6|) . On increasing 
AK (K^ > K^), the microscopic synchronization of 
ensemble 1 is destroyed. 

In summary, we have found new routes to synchrony 
between two AIEOs, characterized by asymmetry in the 
coupling strength. We have also established that the ef- 
fect of white noise is simply to alter the critical values of 
the parameters that control bifurcation, and to change 
the entrainment frequencies, but without otherwise af- 
fecting this route to synchrony. Our main focus has been 
on phase oscillator ensembles, but similar phenomena can 
also be identified in ensembles of other kinds of oscilla- 
tors. 
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FIG. 6: The shift in the entrainment frequencies 0± with 
AK. Here A {1) = A {2 \ B = 0.5, a = 0, Q = and Aw = 1. 
N represents oscillator count in thousands. 
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